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Abstract 

We examine gauge and Yukawa coupling unification in models with gauge-mediated su- 
persymmetry breaking. We work consistently to two-loop order, and include all weak, mes- 
senger and unification-scale threshold corrections. We find that successful unification requires 
unification-scale threshold corrections that are in conflict with the minimal SU(5) model, but 
are consistent with the modified missing doublet SU(5) model for small tan/?, and large tan/3 
with jj, > 0. 
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The apparent unification of the gauge couplings in the minimal supersymmetric 
standard model (MSSM) has sparked much interest in supersymmetric extensions 
to the standard model. In their present form, most phenomenologically viable models 
have two sectors: a hidden sector, in which supersymmetry is broken, and a visible sec- 
tor, which contains the standard-model particles and their supersymmetric partners. 
Supersymmetry breaking is transmitted to the visible sector by gravitational interac- 
tions (as in supergravity-inspired models) or by standard-model gauge interactions (as 
in models with gauge-mediated dynamical supersymmetry breaking). 

Models with gauge-mediated supersymmetry breaking are usually constructed to 
preserve gauge coupling unification to one-loop order. In this letter we will report on 
a closer look at unification in gauge-mediated models. We will present the results of a 
complete two-loop analysis for gauge and Yukawa coupling unification. Our computa- 
tion takes all one-loop thresholds into account, including those at the weak, messenger 
and unification scales. The thresholds include finite terms which turn out to be very 
important for our precision analysis. 

We will present our results in terms of the model-independent unification-scale 
threshold corrections e g and e& ||. These parameters describe conditions that must 
be satisfied by any viable unification model. We will illustrate the range of these 
parameters for the minimal || and (modified) missing-doublet || SU(5) models. 
We will see that present precision measurements exclude the minimal model, but are 
consistent with gauge and Yukawa unification in the modified missing-doublet case. 

In the simplest models of gauge-mediated supersymmetry breaking || , the messen- 
ger sector contains a set of vector-like fields which couple only to a standard-model 
singlet spurion through trilinear terms in the superpotential. The vector-like messenger 
fields are chosen to transform in 5 + 5 or 10 + 10 representations of SU(5). Requiring 
the gauge couplings to remain perturbative restricts attention to at most four 5 + 5 or 
one 10 + 10 plus one 5 + 5 pair of fields. (An additional 5 + 5 pair can be accommodated 
if the messenger particles are sufficiently heavy.) 

We assume that the lowest (S) and highest (F) components of the spurion acquire 
vevs through their interactions with the hidden sector. These interactions remove the 
mass degeneracy of the messenger superfields and transmit supersymmetry breaking 
from the hidden to the visible sector through loop diagrams which contain spurion 
insertions. At the messenger scale, gaugino and soft scalar masses are induced by 
one-loop and two-loop diagrams, respectively. The flavor-blind nature of the gauge 
interactions ensures that flavor- changing neutral currents are suppressed. To this order, 
the soft supersymmetry-breaking A-parameter is not generated. 

The supersymmetric Higgs mass parameter /i and the soft supersymmetry-breaking 
.B-parameter violate a Peccei-Quinn symmetry and cannot be generated by standard- 
model gauge interactions. We will assume that they are generated by some minimal 
mechanism. The region where B = is theoretically appealing |7j because it gives rise 
to a large ratio of vevs (tan (3) without fine tuning. In this region, all CP-violating 
phases are generated only radiatively, so CP violation is naturally small. 

Our approach is as follows. We start with the Fermi constant, Gp, the electromag- 
netic coupling, a em , the Z-boson mass, Mz, the MS strong coupling constant, a s (M z ), 
and the top-, bottom-quark and tau-lepton pole masses, m t , vn^ and m T (for details, see 
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H). We then assume a supersymmetric spectrum and use the full one-loop corrections 
to calculate the DR couplings gi, g%, g$, X t , A& and A r for a given value of tan j3. We run 
these couplings to the messenger scale, M, using the two-loop MSSM renormalization 
group equations. At M we fix the gaugino and soft scalar masses ||. We then run the 
soft parameters back to the squark mass scale, where we impose electroweak symmetry 
breaking and calculate the supersymmetric spectrum. We iterate the procedure several 
times to achieve a consistent solution. 

Our calculations of the one-loop threshold corrections include the finite and loga- 
rithmic terms. The finite corrections, which are often neglected in the literature, allow 
a precise determination of the gauge couplings g% and g<i at the scale Mz @, [K| . The 
finite corrections to the bottom and tau Yukawa couplings also play an important role 
in our analysis. 

Once we determine the gauge and Yukawa couplings at the messenger scale, we 
extrapolate them to the unification scale, Mqut, which we define to be the scale where 
gi and gi meet. We use the usual two-loop beta functions to compute the evolution of 
the gauge and Yukawa couplings. We also include the messenger contributions, those 
listed in Ref. 
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The sum over / runs over all messenger multiplets, n 5 and Uiq are the number of 5 + 5 
and 10 + 10 messenger fields, and 
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At Mqut we set the messenger Yukawas to a common value, y m . We run the messenger 
Yukawas back to the messenger scale according to their one-loop evolution equations, 
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where T = n 5 (3yJ + 2?/|) +n 1Q (6yg + 3y^ + yl) and the Cj's are the quadratic Casimirs, 
3Y 2 /5, 3/4, and 4/3 for fundamental representations. (The messenger- Yukawa evolu- 
tion equations can receive additional model-dependent contributions from the hidden- 
sector particles. The extra terms do not affect the messenger mass splittings, so we 
can ignore them in our analysis. Note that the one-loop equations suffice because the 
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messenger-sector Yukawas enter our calculation only through the messenger threshold 
corrections.) 

From the set of yf(M), we determine the messenger-particle mass spectrum and 
compute the messenger-scale threshold corrections to the gauge couplings, 
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where A = F/S is the supersymmetry-breaking scale and Mf is the messenger fermion 
mass. The second term in the brackets is small for A/Mf <C 1, in which case there is a 
near degeneracy among the masses in the vector-like supermultiplets. Note that there 
are no messenger-scale Yukawa thresholds to this order. 

We iterate this procedure to find a consistent solution in the region between M 
and M GUT . At M GUT we define the threshold corrections for the gauge and Yukawa 
couplings, e g and e&, as follows, 

g 3 (M GUT ) = 9l (M GUT )(l + e g ), 

h(M GUT ) = X T (M GUT )(l + e b ) . (7) 

The parameters e g and e b describe the unification-scale threshold corrections that are 
necessary to achieve unification in any particular model. In what follows, we will 
indicate the allowed ranges of e g and e& for two of the simplest unification models, the 
minimal and the modified missing-doublet SU(5) models. 

In the minimal SU(5) model, the unification-scale gauge threshold correction is 
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where Mh is the mass of the color-triplet Higgs multiplet that mediates nucleon decay. 
Generally, Mh is bounded from below by the proton decay limits [HJ, which imply 
M H >M GUT , so e g >0. 

The missing-doublet model is an alternative SU(5) theory in which the heavy color- 
triplet Higgs particles are split naturally from the light Higgs doublets Q . This requires 
large SU(5) representations, such as the 75 and 50 + 50, so the SU(5) coupling g 5 
diverges below the Planck scale. The modified missing-doublet (MMD) model solves 
this problem for < 1 by lifting the mass of the 50 + 50 to the Planck scale and 
suppressing the nucleon decay rate through an extra Peccei-Quinn symmetry ||. In 
this way the modified missing doublet model can accommodate two color-triplet Higgs 
particles with masses between 10 13 — 10 15 GeV. 

In the modified missing-doublet model, the unification-scale gauge threshold can 
be written as @, ||, [16[ 
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where Mjf is the effective mass that enters the proton decay amplitude, so the previous 
lower bounds on Mh apply here as well. In the MMD case, the effective mass is also 
bounded from above, M|f< 10 20 GeV [§. 
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The Yukawa threshold in minimal SU(5) can be written as follows || [H]], 
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where My is the mass of a superheavy SU(5) gauge boson. For the minimal SU(5) 
model, the most stringent lower limit on My comes from requiring that the 5 + 5 Higgs 



coupling remain perturbative to the Planck scale. This implies My>0.5Mn |[lq| . We 
take the upper limit on My to be the Planck scale, My < 10 19 GeV. 

For the modified missing-doublet model, the Yukawa threshold has the same form 
as eq. (jTU| ), with the color-triplet Higgs mass, M#, replaced by the effective mass, Mjf . 
In this case, the lower limit on My comes from proton decay experiments, which imply 
M v /g GVT >3.8 x 10 15 GeV P[. As before, we impose My < 10 19 GeV. Hence, both 



models have the same upper limit on e&, but the lower limit in the MMD model is 
lower, by virtue of the fact that My can be smaller and Mh larger. 

In what follows, we present our results for gauge-mediated models. In particular, 
we calculate e g , e&, o>gut and M GUT as functions of the input parameters, which we 
take to be tan/5, the numbers and nio, the supersymmetry-breaking scale A, the 
messenger scale M, and the messenger Yukawa at the unification scale, y m . To examine 
bottom-tau unification, we fix the sign of /i to be positive. 

We find the range of olgut and Mqut by scanning over the parameter space, with 
m t = 175 GeV, m b = 4.9 GeV, A < 300 GeV, 1.03 < M/A < 10 4 , 0.03 < y m < 3.0 and 
tan/? in the allowed range. For the case = 1, we determine olgvt — (0.044 — 0.054) 
and M GUT ~ (1.5 - 5.0) x 10 16 GeV. For n 5 = n 10 = 1, we find a GUT ~ (0.062 - 0.28) 
and M GUT ~ (1.2 - 7.0) x 10 16 GeV. 

In Fig. [I] we plot e g and e& for n 5 = 1, M/A = 2, m;, = 4.9 GeV and y m = 1, versus 
A and tan/?, respectively. In (a) we choose tan/3 = 20, while in (b) we take A = 100 
TeV. In each case the short-dashed (long-dashed) lines correspond to a s (M z ) = 0.124 
(0.112). The black bands correspond to a s (Mz) = 0.118 with m t varying from 170 to 
180 GeV. The uncertainty in e& from varying mj = 4.9 ± 0.3 GeV is almost the same 
as that from changing a s (M z ) = 0.118 ± 0.006. 

In Fig. |l](a) we also show the allowed values for e g in the minimal and modified 
missing-doublet SU(5) models. The region of allowed e g in the modified missing- 
doublet model almost completely overlaps the region with a s [Mz) = 0.118 ±0.006. In 
contrast, we see that minimal SU(5) is inconsistent with a s [Mz) by more than 2a. 

For = 1 we find that the messenger sector corrections decrease e g . The change is 
induced by the messenger thresholds and the differences in the two-loop gauge coupling 
evolution. Both of these effects are of approximately equal importance. 

From Fig. [|(a) we see that raising the supersymmetry-breaking scale A decreases the 
size of the gauge-coupling unification-scale threshold. This is because the superpartner 
masses scale with A, and larger masses decrease the size of the required thresholds 
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Figure |Tj(b) illustrates the well-known fact that bottom-tau unification can only be 
achieved for very small (<1.8) or rather large (>35) tan/?. (Very large values of tan/? 
are excluded by the requirement of proper electroweak symmetry breaking.) Figure 
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Figure 1: The unification-scale threshold corrections with n 5 
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(a) The gauge coupling unification-scale threshold 



correction e g versus A, for tan j3 = 20, and a s (Mz) = 0.118 (black band), 0.124 
(short-dashed) and 0.112 (long-dashed), (b) The Yukawa coupling unification- 
scale threshold correction, e&, versus tan /3, for A = 100 TeV and the same val- 
ues for a s (Mz) as in (a). In each case, the black band is obtained by varying 
the top mass from 170 to 180 GeV. The shaded regions indicate the allowed 
range for (a) e g and (b) e& in the minimal and modified missing-doublet SU(5) 
models. 



|T|(b) also shows the allowed bands for e& in the minimal and modified missing-doublet 
SU(5) models. 

As above, we can compare this plot to the case with no messengers. There, one 
typically finds that the bottom and tau Yukawa couplings meet much earlier than the 
scale Mqut, so a rather large and negative threshold correction e& is required. For the 
case at hand, the extra messenger multiplets change the Yukawa evolution equations at 
two loops. More importantly, however, they also increase the gauge couplings, which 
feed into the Yukawa evolution equations and cause the bottom and tau couplings to 
meet even earlier. This makes e& even more negative. 

Fortunately, at large tan f3 there are significant finite threshold corrections to the 
bottom (and to a smaller extent, tau) Yukawa couplings |§. These corrections, which 
are proportional to //tan/3, are sufficiently important to permit bottom-tau unification 
at large tan (3 for // > 0. (The case // < is completely excluded at large tan (3.) These 
finite corrections were omitted in the analysis of Ref. [T2[], which came to a different 
conclusion. 

For ns = 1, the value of e g is not significantly affected by changes in tan /3 or M/A. 
At the smallest values of tan/3, e g increases by about 0.5%, while for M/A = 10 4 , e g 
increases by about 0.2%. The parameter e& is more sensitive to changes in M/A. For 
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Figure 2: The same as Fig. [I], except that n$ = n w = 1, and in (b) A =50 
TeV and M/A = 100. The dot-dashed lines indicate e?, and its lower limit in 
the MMD model for M/A = 2. The line in the upper right-hand corner of (b) 
describes the top-quark threshold, e t , for M/A = 100. 



M/A = 10 4 , the eb curve is 2.5 to 3% higher at intermediate tan/?, and rises to +20% 
at tan (3 ~ 40. 

In Fig. H we plot e g and e b for the case of n 5 = n 10 = 1, versus A and tan/? 
respectively. The other parameters are fixed as in Fig. p], except that in Fig. 0(b), 
A = 50 TeV (to keep the scalar masses unchanged) and M/A = 100. (Two M/A = 2 
curves are shown in dotted lines.) Figure |2|(a) shows that everything shifts because 
of the larger cygut, but the overlap between the band from the MMD model and the 
allowed region for a s [Mz) is still almost complete. In this case, increasing M/A to 10 4 
significantly changes Fig. 2(a). The central value of e g runs from —4% for A = 20 TeV 
to —1.5% for A = 100 TeV. The band for the MMD model is such that the required 
value of e g lies entirely within the band. (Note, however, that n 5 = n w = 1 gives rise 
to nonperturbative couplings above M GUT in the MMD case.) 

The change in Fig. 0(b) as compared to Fig. |I](b) is more dramatic. Because the 
gauge couplings are even larger than in the previous case, bottom-tau unification turns 
out to be barely possible for M/A = 2 (dot-dashed lines). Note, however, that there is 
still a significant region for unification in the missing doublet model with M/A = 100. 
In Fig. ^|(b) we also show the necessary threshold, et, for top-tau Yukawa unification. 
We see from Fig. |2|(b) that the top, bottom and tau couplings unify at the largest 
values of tan/3 (in the region where B ~ 0). Such a unification is expected in SO(10) 
models. However, the thresholds in any particular SO(10) model must be calculated 
to be sure the model is consistent with data. 
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